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We study the finite temperature, low energy, long wave-length spectrum of the dynamic structure 
factor of the spin-1/2 antiferromagnetic Heisenberg chain in the presence of exchange anisotropy 
and external magnetic fields. Using imaginary-time quantum Monte-Carlo we extract parameters, 
relevant to characterize a renormalized Luttinger liquid. For small momentum our results are con¬ 
sistent with a change from propagating spinon density waves to spin diffusion, described by a finite- 
frequency spin-current relaxation rate. Results for this relaxation rate as well as other Luttinger 
liquid parameters are presented versus temperature, momentum, magnetic field, and anisotropy, in¬ 
cluding finite-size analysis, and checks for anomalous diffusion. Our results are consistent with exact 
diagonalization and Bethe Ansatz, where available, and with corroborate findings of other previous 
studies using bosonization, transfer matrix renormalization group, and quantum Monte-Carlo. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.50.Ee, 75.40.Mg 


I. INTRODUCTION 

The one-dimensional (ID) spin-1/2 antiferromagnetic 
Heisenberg chain (AFHC) is one of the most fundamen¬ 
tal models of quantum many-body physics. It is rel¬ 
evant to low-dimensional magnets^i ultra-cold atoms?^ 
nanostructures)^ and - seemingly unrelated - fields such 
as string theory^ and quantum Hall systems^ In the pres¬ 
ence of an external magnetic field, its generalization to 
anisotropic exchange, the XXZ model, reads 

H = jY.[^{S+S-^,+S-St^,)+^SfSt^,-hSt] , (l.I) 

where J > 0 is the antiferromagnetic exchange interac¬ 
tion with an anisotropy ratio A, are the spin opera¬ 
tors on site I of a chain of length L with periodic bound¬ 
ary conditions (PBC), and B = g^Lshh is the magnetic 
field. 

Experimentally, dynamical correlation functions of the 
AFHC have recently become accessible to a variety of 
high-resolution spectroscopies at finite temperature and 
in the presence of external magnetic fields, e.g. inelastic 
neutron scattering (INS))^^— high-field nuclear magnetic 
resonance (NMR))S“— muon spin-resonance (pSR))^^ 
and magnetic transport 

Theoretically, and while the AFHC is integrable, its 
dynamical spin correlation functions remain a major 
challenge. Analytically, signihcant progress has been 
made at low temperatures by calculating multi-spinon 
response functionsii from Bethe Ansatz, including the 
cases of h 7 ^ 0 and A 7 ^ 1. Important insight has also 
been obtained in the continuum limit and at low tem¬ 
peratures by bosonizatioiii ^^d^ Perturbation theory al¬ 
lows to access regimes of A ^ li^ Numerically, at h- 
nite temperature, time-dependent density-matrix renor¬ 
malization group (t-TDMRG)^ is a very powerful ap¬ 
proach. However, at present, entanglement growth re¬ 
mains a limiting factor on the time window to access 
low-energy, long-wavelength dynamicsj ^^i^^ Recently, dy¬ 


namical quantum typicality (t-QT) has been shown to 
overcome such limitations)^^— however, only at high 
temperatures. Early on, exact diagonalization22“— and, 
more recently, Lanczos variants2£ have been used. How¬ 
ever, finite-size or Krylov-space-dimension effects limit 
their spectral resolution. 

Quantum Monte-Carlo (QMC) is an additional com¬ 
plementary approach for the HAFC. It is applicable from 
low to high temperatures and allows to consider systems 
almost in the thermodynamic limit, including finite mag¬ 
netic fields and anisotropy. Static correlation functions 
can be obtained from it with arbitrary precision, however, 
the evaluation of spin spectra from QMG^“— requires an¬ 
alytic continuation of imaginary-time data, leading to er¬ 
rors from maximum-entropy approaches— This renders 
the evaluation of the dynamic structure factor S{q, lo) by 
QMC very challenging, in particular for small momenta 
g <C 1 , where sharp low-energy spectral features are ex¬ 
pected. Recently therefore, QMC results^^ for S'(g, r) at 
imaginary time r have been compared directly to sug¬ 
gestions from bosonization and time-dependent transfer- 
matrix renormalization group (t-TMRC) )^^i^^ corrobo¬ 
rating a picture of finite-temperature diffusive spin dy¬ 
namics at finite frequencies in the low-energy, long- 
wavelength limit. Such spin diffusion in the HAFC is a 
long-standing issuO)^ relating the magnetization dynam¬ 
ics to the question of dissipation of spin currents and the 
quest for a spin-Drude weight in the XXZ model, see e.g. 
Ref. [ 2 ^ and references therein. 

In this context, the aim of the present work is to extend 
the information obtained from the QMC approach of Ref. 
1 ^ into various directions, including physical as well as 
technical aspects. The manuscript is organized as follows. 
In Sec. m we describe our method of analysis. In Sec. IHII 
we revisit Ref. from three directions, first, by subject¬ 
ing its findings to a finite-size scaling analysis, second, 
by comparing our extended results with exact diagonal¬ 
ization, and third, by checking for anomalous corrections 
to diffusion. The following Secs. IIVI and |V] extend the 
analysis in two additional ways, namely, by studying the 
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impact of finite magnetic fields and anisotropy. Section 
ED summarizes our findings. 


II. QMC APPROACH 

The prime object of interest in this work is the Fourier 
transform of the longitudinal dynamic structure factor 


S{q,U}) 




( 2 . 1 ) 


at small momentum q and frequency w, and its 
corresponding retarded dynamical spin susceptibility 
Xret( 9 ,w) = x' + *X" is related to S{q,uj) by 
the fluctuation-dissipation theorem x"(( 7 ,w) = [1 — 
exp(—/3w)]S'(g,w)/2 at inverse temperature /3 = l/T. 
Here t refers to real time and S{q,t) = with 

Sg = being the spin z component at momen¬ 

tum q. 

S{q,u;) is related to the imaginary-time structure 
factor S{q,T) = {Sy{T)Stq) by analytic continua¬ 
tion through the usual integral transform S'(q,r) = 
dujK{uj,T)S{q,u})/TT with a kernel K{u),t) = e~'^^ + 
g-(/3-T)i^_ While QMC allows to evaluate S{q,T) with 
high precision, inverting the latter integral transform, 
e.g. by maximum entropy methods,— is mathematically 
ill-posed and introduces errors which can be consider¬ 
able in particular at small q and w— A central point 
of the present work is, that instead of performing such 
analytic continuation, we follow^^ and express x{q:<^) 
through only a few relevant parameters, by fitting nu¬ 
merical imaginary-time QMC data to an educated ana¬ 
lytic guess for r). While this avoids errors from an¬ 
alytic continuation, it obviously represents a bias which 
requires posteriori quality checks. 

Throughout this paper we will be interested in the 
long-wavelength, low-energy limit, and in a tempera¬ 
ture range where we may assume that the physics of Eq. 
(11.11) can be described by a renormalized Luttinger liquid 
(LLQ). Therefore we start from 


Xret (g, W) 


Kvq^/2 tt 

^2 _ y2q2 _ w) ■ 


( 2 . 2 ) 


Except for Hret, this form of Xret is dictated by the re¬ 
sponse of a free LLQ with Luttinger parameter K and 
spinon velocity v. All deviations from the free LLQ 
are encapsulated in the irreducible density self-energy 
nret( 9 ,w). Causality requires that the real(imaginary) 
part of the latter is an even(odd) function in w. Assum¬ 
ing that nret( 9 ,w) is non-singular for small w, and since 
we are interested in the low-w limit, we expand to O(w^) 

nret(g,w) ~ iio(g) +iui{q)u} -f U 2 [q)uj‘^ + ... (2.3) 


with Ui{q) G M. Expanding also to O(q^) for g <C 1, 
inversion symmetry forces 

Hret)?, w) ~ cv'^q^ — 2ixoJ — hoj'^ -I- ... (2.4) 


with real constants c, 7 , and 6 , and all terms of order 
higher than uP, quj, and q^ have been dropped. Since 
we will confine ourselves to |A| < 1 and the longitudinal 
dynamical structure below the saturation field h < he = 
2 we may discard terms constant in q and ui in Eq. (12.411 . 
since they would lead to a gap in Xret( 0 ,a;). 

We emphasize, that Eq. p.4ll follows solely from 
causality, analyticity, and symmetry. Therefore, bosoniz- 
ing the HAFC and performing proper perturbation the¬ 
ory beyond the free LLQ must result in exactly this small 
q, u j expansion for H^et ( 9 , w). This is consistent with Refs. 
[l8l.fl^l3^ an din where explicit expressions for c, 7 , and 
b have been obtained to second order in Umklapp scat¬ 
tering and first order in band curvature. In contrast to 
this perturbation theory however, the spirit of our work 
is very different. In fact we will view v, K, c, 7 , and b as 
parameters to be determined by fitting to imaginary-time 
QMC results for x( 9 ,t). This is a non-perturbative ap¬ 
proach. To this end we map Eqs. (12.2L (12.41) onto bosonic 
Matsubara frequencies ujn = 2'KnT with integer n. 


x{q,^n) 


_ KqVqq^/{2TT) _ 

(1 -k b)ujl -k (1 -k c)u2g2 + 27 ,|a;„| ’ 


(2.5) 


and transform this to imaginary time with 


X ( 9 , 't) = 2 ^ cos(a;„T)x(9, w„) - x(9, 0 ) , ( 2 . 6 ) 

n—0 

Several comments are in order. First, the static sus¬ 
ceptibility Xq = Xi?(9i0) of the XXZ model is known 
to vary with q, e.g., at A = 1 and h = 0 it increases 
monotonously as 9 —>■ tt/2. However, Xq = K/[2ttv(\ -k 
c)] resulting from Eqs. (j2.2l) . (12.41) is momentum- 

independent. To fix this shortcoming of the free LLQ 
theory, we allow for an additional momentum dependence 
K Kq, V ^ Vq in Eq. (12.51) - albeit weak at 9 <C 1 -, 
when matching Eqs. (I2.2L (12.41) with QMC. Second, the 
constants b and c in Eq. (j2.4l) are redundant when fit¬ 
ting to QMC because they can always be absorbed into 
a renormalization of Kq, Vq, and 7 . Without loss of gen¬ 
erality we therefore fix b and oM- to the values obtained 
from the perturbation theory of Refs. [llElsD, andS. 
Third, as an additional generalization of Eq. (12.4|) we al¬ 
low for a momentum dependence 7 —>■ 7 ^. Fourth, Eq. 
(12.21) with Eq. (12.4|) does not capture the finite width of 
the spectral function x"( 9 , w) at T = 0, resulting primar¬ 
ily from the two-spinon continuum. However, at q/i: <C 1 
the latter width is of order Jq^, which is negligible against 
Xq for those wave vectors and temperatures which we will 
be interested in. 

In summary, our procedure consists in evaluating 
x( 9 ,t) by QMC and fitting the result at each q, using 
Eq. (j2.6p . in terms of three numbers: Kq, Vq, and Xq- 


A. Spin Diffusion 

For hnite Xq to 7 ^, Eq. (12.2p displays a 
diffusion pole with a diffusion constant Fg = (1 -k 
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FIG. 1. (color online) Differences between the fitted curves 
and the original QMC results at the isotropic point for differ¬ 
ent wave vectors and temperatures. The system-size consid¬ 
ered here is L = 256. Error bars shown in the figure indicate 
the statistical noise of the QMC data. 


c)v^/{2^q). Therefore the prime effect of the renormal¬ 
ization ni.et( 5 , w), which we study in this work, is to gen¬ 
erate a hydrodynamic regime in which the magnetization 
dynamics of the HAFC is diffusive. This also has conse¬ 
quences for the regular spin conductivity areg{<^), which 
we sketch briefly next. 

The spin current is defined by the lattice version of the 
continuity equation dtS^ = [1 —exp(—ig)] and its con¬ 
ductivity a(uj) = (l—e^‘^)Re can 

be decomposed as 

(t(w) = T>(5(w) -I-(Treg(w) , (2.7) 

where D is the enigmatic spin Drude weight and (Treg(w) 
is the regular part of the spin conductivity. The Drude 
weight of the HAFC has been of interest for more than 
two decades by now. We refer to Ref. [2^ for the recent 
status. The regular conductivity is related to the spec¬ 
trum y"(g, w) trough the continuity equation 

crreg(w) = liin -^x"(g,w). (2.8) 

If x(g, cu) exhibits a diffusion pole in the frequency do¬ 
main, then in the time domain Eq. 12.81 implies exponen¬ 
tial decay of the current on a time scale I /70 given by 
the current relaxation rate 70 In the frequency domain 
it implies a Lorentzian line shape of (7reg(uj) with width 
70 - We emphasize that Eq. (j2.8ll provides no insight into 
the value of D. Coexistence of a diffusive regular con¬ 
ductivity with a finite Drude weight has been suggested 
based on bosonization, t-TMRG, and memory-function 
methodsi^i^ But also strong momentum dependence be¬ 
yond Eq. (j2.4p has been invoked to question this, based 
however on small systems at high temperaturesi^ For the 
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FIG. 2. (color online) Finite-frequency spin-current relax¬ 
ation rate jg of the isotropic Heisenberg spin chain versus 
inverse system-size 1/L and momentum q for chain lengths 
64 < L < 256 at two temperatures T/J = 0.1, 0.25. For mo¬ 
mentum dependence of 7 ,, the largest system-size L = 256 is 
chosen. 


remainder of this work we will refrain from speculations 
on D. 


III. ISOTROPIC CHAIN 

In this section we extend previous work^ based on the 
approach described in Sec. im into several new directions. 
First we study the effects of system-size. Then we per¬ 
form a consistency check of our approach by comparing 
to available results from exact diagonalization. Finally, 
we test if the current relaxation rate shows relevant en¬ 
ergy dependence. 


A. System-size and momentum dependence 

In Ref. [ 3 ^ the approach described in the previous sec¬ 
tion has been applied to the isotropic HAFC, i.e. at A = 1 
and h = 0. That work has proven feasibility of the 
method. In particular it was shown, that Kq and vq 
extracted from it agree very well with thermodynamic 
Bethe Ansatzj ^^d^ as well as that the relaxation rate 70 
obtained was consistent with that from bosonization and 
t-TMRG— to within factors of 2. We will not repeat 
this analysis here. Instead we note that the findings of 
Ref. m where based on a single system-size, i.e. L = 128, 
and it remained unclear if parts of its results, e.g. the 
momentum dependencies observed for 7 ^ were hnite-size 
effects. Therefore, in this subsection, we will analyze 7 ^ 
versus system-size. 

To begin we comment on the accuracy of the fitting 
procedure summarized in the last paragraph before Sec. 
Ill Al To this end Fig. [T] shows the minimum difference 
between our QMC results and Eq. (12.61) for optimized 
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FIG. 3. (color online) Hydrodynamic spectra in the frequency 
domain with parameters fitted to QMC as compared to ED, 
both for length L = 18, momentum q = tt/Q, the isotropic 
point A = 1, and three temperatures TlJ = {0.33,0.5, 1.0}. 


parameters itTg, fg, and 7 g, for different temperatures 
and small, nonzero wave vectors on a system with 256 
sites. In addition each panel also shows this difference 
assuming that 7 g = 0, i.e. discarding spin diffusion. The 
figure also displays the statistical noise of the QMC data 
as error bars. Several points have to be mentioned. First, 
the quality of our imaginary-time QMC results is worth 
noting by observing the absolute order of magnitude on 
the y axes of this figure. Second, for all panels shown, 
including a finite ■jq leads to clearly better agreement 
with QMC than setting 7 g = 0. We emphasize, that 7 
is ^ 0( J/20... 100), which implies that our method is 
capable to sense spectral structures on such small energy 
scales. This is unlikely for standard maximum-entropy 
analytic continuation. Third, at intermediate tempera¬ 
tures, T = 0.25 J, using Eqs. ( 1 ^ . UM with finite 7 ^ is 
well within the standard deviation of QMC for all mo¬ 
menta depicted. For lower temperatures this is so at 
least for the smallest momenta, i.e. as q —>■ 0. There¬ 
fore these plots also show that, while a finite 7 ^ clearly 
improves agreement with QMC, the range of validity in 
momentum space of the purely hydrodynamic, diffusive 
description shrinks as the temperature is lowered. Simi¬ 
lar observations have been made in Ref. IdO . Future work 
should elaborate on this. Finally, the figure shows that 
7 g decreases with temperature 

Next we turn to the system-size dependence of 7 ^. 
To that end we consider four chains with sizes of L = 
64,128,192, 256. The results are shown in Fig. [5] a), b) 
for two different temperatures T/J = 0.1, 0.25 and two 
momenta. Note that these momenta, q = 7r/32, 37r/32, 
have been chosen to be identical despite the differing 
system-sizes, which implies that they are not necessar¬ 
ily the smallest wave vector possible for a given L. The 
main point of this figure is, that we observe no finite- 
size dependence. This is even more so in view of the y 
axes scales depicted. Therefore our results can safely be 
considered as being in the thermodynamic limit. 

In panels c), d) of this figure we show 7 ^ versus mo¬ 
mentum for the largest system, i.e., L = 256. Obviously, 
for the lower temperature T = 0.1, 7 g slightly decreases 
with This is fully identical to the findings of Ref. 


FIG. 4. (color online) Exponent a of the spin-density relax¬ 
ation rate 701 “ of the isotropic Heisenberg chain versus a) 
inverse system-size 1/L at momentum q = 7r/32 of system- 
sizes 64 < L < 256, b) momentum q for a system-size of 
L = 256. In both panels, two temperatures TlJ = 0.25, 0.1 
are considered. 

[ 3 ^, where however L = 128 was used. I.e. this decrease 
is not a finite-size effect as was speculated in the latter 
work. Rather, and in view of the deviation in Fig. [T]e), 
f), we believe that the crossover from a constant to a 
weakly momentum-dependent 7 ^ is another indication of 
a shrinking hydrodynamic regime as the temperature is 
lowered. 


B. Comparison with exact diagonalization 

Apart from Fig.[Tl another quality check of our method 
can be constructed by comparison with results from exact 
diagonalization (ED). While the latter are beyond any 
doubt, they can only be obtained on small systems, in 
particular if complete diagonalization is performed using 
a canonical ensemble. Small system-sizes imply that the 
smallest, nonzero wave vectors remains rather large, and 
moreover, that acceptably smooth spectra can only be 
obtained for rather large temperatures. Nevertheless, it 
is very tempting to compare our QMC approach with ED 
results. 

Figure [3] compares the evolution with temperature of 
the longitudinal structure factor S{q,uj)^ obtained from 
both, the QMC method of Sec. HIl and ED on chains of 
identical length, i.e. L = 18, at the smallest, nonzero 
wave vector q = tt/9. This comparison is remarkable. 
First, it is obvious, that as T is lowered, the agree¬ 
ment between ED and the QMC approach improves 
monotonously. At the lowest temperature the width of 
the two spectra are almost identical. The peak heights 
are still different. We attribute this to two effects. First, 
at low temperature respective finite-size effects are more 
prominent. Second, the width of the ED spectrum is still 
slightly larger than that of QMC. In turn, sum-rule effects 
will decrease the ED’s peak intensity. As we increase T 
the line shape of the ED clearly broadens beyond that of 
the QMC aXTlJ = 0.5, and for T/J = 1.0 the agreement 
is lost. 

Exactly this variation with temperature is expected. 
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FIG. 5. (color online) a) Spin current relaxation rate 7 q, 
b) Luttinger parameter K, and c) spin velocity v at finite 
temperatures T/J = 0.1, 0.25 obtained from QMC of a chain 
with 64 sites. Luttinger parameter K and spin velocity v 
at zero temperature known from BA in terms of magnetic 
field are also shown. In the case of QMC, the momentum is 
q = 'kI32. 

For low T a description of the HAFC in terms of a renor¬ 
malized LLQ shonld hold, explaining the rather similar 
spectra from ED and our QMC approach. We speculate 
that this agreement should even improve for T < 0.33 J, 
where however ED spectra are too noisy. Increasing the 
temperature, it is known that up to T/J < 0.1 per- 
turbatively renormalized LLQs provide a good descrip¬ 
tion of thermodynamic properties of the HAFC which 
agree rather well with static QMC or TBAfii^ while for 
T/J'> 0.25 the agreement starts to deteriorate. In view 
of this breakdown of the LLQ description, it is natu¬ 
ral that ED and our QMC approach start to differ for 
T/J > 0.5, as in Fig. [31 



1ir/32 3jr/32 5jr/32 7ir/32 

q 

FIG. 6. (color online) Momentum dependence of the spin 
current relaxation rate for two temperatures a,)T/J = 0.25, 
h)T/J = 0.1, and different external magnetic fields is plotted. 
The system-size considered here is L — 64. 

This claim is consistent with numerical simulations 
of spin-density dynamics on the basis of classical 
mechanics 1 ^ 

Allowing for the modification of Eq. (|3T]) within our 
approach is straightforward and leads to one additional 
fitting parameter. The results for a at two temperatures 
T = O.IJ, 0.25 J are summarized in Fig. [4]. Panel a) of 
this figure depicts finite-size scaling of a for system-sizes 
64 < L < 256 and one momentum q = 'r:lZ2. a versus 
momentum for a system-size of L = 256 is shown in panel 
b) of this figure. To summarize, a. remains very close to 
one, i.e. any tendencies towards super(sub)diffusion at 
low temperatures are negligible and can safely be dis¬ 
carded. We note in passing, that all modifications we 
observe in Kq, Vq and 7 g, due to allowing for a ^ 1, are 
also negligible. 


C. Anomalous diffusion 

As a final check we test the stability of our approach 
against modifications of the basic diffusion law by su- 
per(sub)diffusion. The latter is defined by a density 
which spreads in time according to (r^) ~ with 
a > 1(< 1). In terms of Eq. (12.41) . this implies a sub¬ 
stitution 

7 ^ —>• 7w“ (3.1) 

with a 1. Recently, anomalous diffusion has been 
claimed to result from Lindblad quantum master equa¬ 
tions at the isotropic point and infinite temperatures 4^1^ 


IV. MAGNETIC FIELD DEPENDENCE 

In this section, we will extend our approach to the 
case of finite magnetic fields in the isotropic chain. Since 
we consider the longitudinal dynamical structure factor, 
the spectrum remains gapless below the saturation field 
he = 2 . 

First, Fig. |S^) shows jq versus field for two temper¬ 
atures T = 0.1, 0.25 J at wave vector q = 7r/32 in a 
system with L = 64. This evidences a monotonous de¬ 
crease of the relaxation rate with magnetic field. Such 
a behavior is consistent with predictions at low fields 
from bosonization4^ While overlap between the spin cur¬ 
rent and the conserved heat current of the HAFC at 
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FIG. 7. (color online) Differences between the fitted curves 
and original QMC data at momentum q = 7r/32 and two 
temperatures TjJ = 0.25, 0.1. In panels a) and b), magnetic 
field is set to h = 0.25 and anisotropy is set to A = 1 while 
in panels c) and d), h = 0 and A = 0.8. The system-size for 
all cases is L = 64 and the error bars indicate the statistical 
noise of the QMC data. 

nonzero h irrevocably imply a finite spin Drude weight for 
0<h<hc^^— this allows for no conclusion on the regular 
conductivity (Ti.eg(w) from Eqs. (12.71) . (|2.8D . In Refs. 
and m this has been dubbed coexistence of ballistic and 
diffusive transport channels, meaning that in the time 
domain, and on relatively short scales, currents relax 
diffusively onto a constant infinite-time limiting value. 
Therefore, while our nonzero "fq does not contradict a fi¬ 
nite Drude weight, it is nevertheless reassuring that the 
decrease of we observe is indicative of a ‘more’ bal¬ 
listic behavior of the system as h increases. 

Next we turn to the remaining two fitting parameters 
Kq and Vq. As the magnetic field sets the chemical poten¬ 
tial of the Jordan-Wigner fermions, it is clear that apart 
from changes in jq versus h we also expect the LLQ pa¬ 
rameters Kq and Vq to be sensitive to h. This is shown 
in Fig. Oa) and c) for the same wave vector and temper¬ 
atures as in panel a). At T = 0 the LLQ parameters 
can be obtained from Bethe Ansatz (BA)<^ Regarding 
our approach it is very satisfying to see that Kq and Vq 
almost quantitatively agree with BA, and moreover that 
the remaining discrepancies diminish as T —>■ 0 in the 
QMC approach. We emphasize that it should come as 
no surprise, that expectation values obtained from static 
QMC agree with other exact methods. Rather the main 
point is that Kq and Vq within our QMC approach follow 
from hts to a dynamic structure factor. We note that 
the downturn of Vq with h is consistent with the phys¬ 
ical picture of a cosine band being filled completely as 
h —y he- 

Figure [HI displays jq versus momentum for two temper¬ 
atures T = 0.1, 0.25 J and various magnetic fields h = 0, 
0.5, 1.0, up to half of the saturation field. This figure 
shows, that increasing h decreases the slope of 7 ^ and 
modifies the curvature as well. However, ^q(h > 0) is not 


FIG. 8. (color online) Anisotropy dependence of the spin- 
current relaxation rate obtained from bosonization (solid 
curves) and our method (symbols) for two temperatures 
T/J = 0.1, 0.25. For QMG, the first nonzero momentum 
g = 7r/32 on a system of length L = 64 is considered. 


less then 7 q(h = 0) for large momenta. To gauge if the 
introduction of 7 g improves the agreement between the 
renormalized LLQ and the imaginary-time QMC data, 
as compared to the free LLQ, also for finite magnetic 
fields, we display the difference between x(q, t) from Fq. 
p.5l) and that obtained from QMC in Fig. [Tji) and b). 
This figure clearly supports 7 ^ ^ 0. Similar however to 
the zero-field case, the QMC data also suggests that at 
low temperatures, additional, albeit smaller, renormal¬ 
izations beyond diffusion are present. 


V. ANISOTROPY 

Finally, we consider the role of finite anisotropy for 
the range of 0 < A < 1. To this end Fig. |8] dis¬ 
plays the relaxation rate versus A for two temperatures 
T/J = 0.1 and 0.25 obtained from two methods, i.e. our 
QMC approach and the perturbation theory based on 
bosonization i^^'^^ For QMC the results are plotted for 
the smallest, nonzero wave vector on a chain of length 
L = 64, i.e. q = 7r/32. The bosonization result is of 
the form 7 (T) = /(AT, where / is a function of 

the Luttinger parameter K and the spin velocity v which 
are exactly known from BA at zero temperature. This 
figure shows, that the QMC results behave qualitatively 
very similar to the findings from bosonization yet, QMC 
is consistently larger than the latter. The figure does 
not show the point A = 1, since at that point a similar 
comparison has already been performed in Ref. 1^ . The 
evolution of 7 ^ with A as the latter approaches zero is 
very much consistent with the fact that at A = 0 the 
HAFC is in the XY limit, i.e. a free Fermi gas, which 
displays fully ballistic spin transport, i.e. 7 = 0 . 

As for the case of finite magnetic fields, in Fig. [S] we 
now turn to the dependence of the two LLQ parameters 
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FIG. 9. (color online) a) Luttinger parameter and b) spin 
velocity at finite temperature obtained from QMC in terms 
of anisotropy for the first nonzero momentum at chain length 
L = 64. Luttinger parameter and spin velocity at zero tem¬ 
perature known from BA are also shown for comparison pur¬ 
poses. 


Kq and Vq for the smallest wave vector on a L = 64 site 
chain. Once again, values for these parameters at T = 0 
are known from BA, i.e. 


K = -—, V = - -— . (o.ij 

TT — arccos A 2 arccos A 

In Fig. [5] results from the fits to QMC are compared with 
BA. As for the finite field we observe very good agree¬ 
ment, which improves as T —>■ 0. 

As for the previous sections, we also consider the ef¬ 
fects of anisotropy on the momentum dependence of the 
relaxation rate. This is shown in Fig. [TU] at two tem¬ 
peratures T = 0.1, 0.25 J for the isotropic point A = 1 
and two values of anisotropy A = 0.2, 0.6. Similar to the 
finite-field case, the slope and curvature of 75 (A) versus q 
varies with A, albeit less strong. For the higher tempera¬ 
ture, i.e. T = 0.25 J, only a weak momentum dependence 
can be observed at the isotropic point and an increasing 
function in the anisotropic case. For T = 0.1 J, 75 de¬ 
creases with g at A = 1 , while for the anisotropic case, 
this behavior is reversed. 

The quality of the QMC fits to Eq. (12.61) for a finite 
A = 0.8 are depicted in Fig. [T] They exhibit the same 
trend as in all similar comparisons discussed in this work, 
i.e. a finite 7 clearly improves agreement with QMC as 
compared to 7 = 0 , with however systematic deviations 
still present for very low temperatures and hnite mo¬ 
menta. 


FIG. 10. (color online) Momentum dependence of the spin- 
current relaxation rate 7 ^ for different anisotropies A = 
1.0, 0.6, 0.2 at two temperatures a) T — 0.25J and b) T = 
O.IJ. The system-size here is set to L = 64. 


VI. SUMMARY 

In summary, we have explored the dynamic struc¬ 
ture factor of the antiferromagnetic Heisenberg chain at 
long wave lengths and finite temperatures. For this we 
have used the response function of a Luttinger liquid in¬ 
cluding self-energy corrections which were parametrized 
in terms of only a few relevant parameters. We have 
then determined these parameters by fitting to high- 
precision imaginary-time quantum Monte-Carlo calcula¬ 
tions. Quite generically we find, that the self-energy 
renormalizes the low-energy, long wavelength density 
spectrum of the Luttinger liquid from a propagating 
‘sound’ mode into a diffusion mode. Apart from the 
standard Luttinger parameters, the diffusion mode in¬ 
troduces one additional parameter, i.e. a current relax¬ 
ation rate. We have analyzed this relaxation rate with re¬ 
spect to several physical variables: temperature, momen¬ 
tum, magnetic field, and anisotropy. Moreover, we have 
shown consistency of our approach by performing hnite- 
size analysis, checking for anomalous diffusion, and com¬ 
paring with exact diagonalization and results from Bethe 
Ansatz. Our findings are qualitatively consistent with 
previous studies using perturbation theory and bosoniza- 
tion, as well as with a limited QMC analysis, considering 
only the role of temperature and momentum for a single 
system-size. Future work should enlarge the variational 
space for the self-energy, in order to capture corrections 
to purely diffusive behavior. 

Note added: After completion of this work we have be¬ 
come aware of very recent results^ obtained in a differ- 

















ent context, with different methods, which are also con¬ 
sistent with a finite zero-frequency limit of the regular 
spin conductivity in the HAFC. 
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